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rr' '. Abstract 
> 

^O ■ We present an experimental and numerical study of the effects 

^^ I of decoherence on a quantum system whose classical analogue has 

i^ ■ Kolmogorov-Arnol'd-Moser (KAM) tori in its phase space. Atoms are 

f^*) I prepared in a caesium magneto-optical trap at temperatures and den- 

Q^ ' sities which necessitate a quantum description. This real quantum 

0^ I system is coupled to the environment via spontaneous emission. The 

r~| ' degree of coupling is varied and the effects of this coupling on the 

fS il quantum coherence of the system are studied. When the classical dif- 
fusion through a partially broken torus is < h, diffusion of quantum 

Ch I particles is inhibited. We find that increasing decoherence via sponta- 

2 ' neous emission increases the transport of quantum particles through 

^^1 the boundary. 



1 Introduction 

The study of decoherence in a quantum system has been a subject of much 
interest in recent years. Since the emergence of quantum mechanics some 
seventy years ago, a central problem in its interpretation has been the fact 
that the linearity of the Schrodinger equation allows macroscopic physical 
states to be derived from arbitrary superpositions of other states, a situation 
which is quite clearly in contradiction to our experience. This was originally 
illustrated in the familiar paradox of Schrodinger's cat. One simple explana- 
tion is to suppose that an initially classical state for a macroscopic system 
will always evolve into another classical state so that for a universe with 
a suitable initial condition, non-classical states occur only when carefully 
prepared by some experimentalist. In this interpretation Schrodinger's cat 



presents no paradox, but merely lies outside our normal realm of experi- 
ence. This hypothesis can be rejected in the light of quantum analysis of 
real macroscopic, classically chaotic systems. Calculations (see for example 
Zurek's discussion in terms of celestial objects |1|) indicate that in these 
systems, quantum dynamics differ from classical after relatively short times, 
and lead to flagrantly non-classical states. Any successful theory must ex- 
plain why these states are not found in practice. 

There have been several more sophisticated explanations as to why quan- 
tum states are not evident in our everyday world. The first was Neils Bohr's 
Copenhagen Interpretation (CI) [Q|. Bohr theorized that quantum theory 
was not universal — a classical apparatus was necessary to carry out the 
measurements. He also required that the boundary between quantum and 
classical domains was mobile and could be pushed back by appropriate ap- 
paratus. Another significant school of thought was Hugh Everett's Many 
World's Interpretation (MWI) [§,§]. He proposed that each time an inter- 
action which would produce a superposition takes place, the wave-function 
of the universe bifurcates to give an ever increasing number of 'branches'. In 
1952, David Bohm ||5| detailed an interpretation based on 'hidden' variables 
which had previously been investigated to some extent by de Broglie [Hj. 
Bohm theorized that these hidden variables would allow "a detailed causal 
and continuous description of all processes, and not require us to forego the 
possibility of conceiving the quantum level in precise terms." Many other 
approaches to the problem have been taken, examples of which can be found 
in&i. 

A more recent and arguably more satisfactory solution to the problem of 
quantum-classical correspondence has been championed by Wojciech Zurek 



and co-workers [0, 11]. They suggest that any real, open quantum sys- 
tem leaks coherence to its surroundings via extraneous degrees of freedom 
which are coupled to the environment. A suitable environment might be 
provided by air molecules, cosmic background microwave photons, or the 
vacuum point fluctuations of the electromagnetic field. The rate at which 
this decoupling proceeds depends on the particular state of the system, the 
dynamics of the system and the form of the interaction with the environ- 
ment. The states most resistant to decoherence form a (not necessarily 
orthogonal) basis of 'pointer' states which are 'selected'; the density ma- 
trix rapidly becomes diagonal in this basis |l^. Any particular realization 
of the system quickly evolves into one of these basis states, and only mea- 
surements in this basis can be thought of as revealing the true state of the 
system. It appears to be a general result that the "pointer" states are the 
"most classical" states of the system. In most cases for macroscopic systems 



the pointer states are expected to be those which are locahzed to a high 
degree. This arises from the fact that in most cases the system's interac- 
tions with the environment are chiefly position dependent. Superpositions 
of the pointer states of a system will collapse on a time scale determined by 
the coupling. In a measurement of a microscopic system by a macroscopic 
observer, the wave-function collapse of traditional measurement theory is 
caused by decoherence of the meter entangled with the system. Bizarre 
macroscopic quantum states might in theory be prepared, but will survive 
only for a vanishingly short time, and the classical description of the world 
as we observe it will be regained. Quantum mechanics without fundamental 
modification retains its position as the true description of the universe, and 
the limited set of states which we observe around us is accounted for. 

In this paper, we study experimentally the effects of increasing the cou- 
pling of a quantum system to its environment. Our experimental system is a 
low density cloud of super cool (~ 15fiK) caesium atoms prepared in the ul- 
tra high vacuum glass cell of a magneto-optical trap (MOT) . Because of the 
low densities achieved in a MOT, the interactions of the atoms are negligible 
and each atom can be modeled as a quantum particle in a periodic poten- 
tial (see Section H). With 10^ atoms trapped per experimental run, we deal 
with statistically significant numbers of particles so the experimental dis- 
tributions will closely approximate quantum probability distributions. We 
temporally modulate a standing wave optical potential which then creates 
Kolmogorov-Arnol'd-Moser (KAM) tori (impenetrable momentum barriers) 
in the classical phase space of the system. For an increased Rabi frequency 
(i.e. increased perturbation), holes appear in the barriers, now called can- 
tori, through which the atoms can diffuse. The diffusion rates of classical 
and quantum particles are distinctly different, with quantum diffusion be- 
ing largely suppressed by the cantori. The decohering effects of coupling to 
the environment bring the behaviour of the quantum ensemble towards the 
classical limit. 

2 The Kolmogorov-Arnol'd-Moser Theorem 

For a closed system with an integrable Hamiltonian, all solutions lie on two 
dimensional tori embedded in 2N — 1 dimensional phase space, where A^ is 
the number of degrees of freedom of the system. Each solution or trajectory 
is indefinitely confined to its own torus. The motion of a trajectory on a 
torus can be described by two coordinates, 6i and 62 on [0, 2it). The winding 



number w of the torus is defined as 

w = — , where Uj = 9j, for j = 1, 2. (1) 

If the two frequencies ujj are commensurate then the winding number of the 
torus is rational, and trajectories on the torus are periodic. If the frequencies 
are incommensurate the winding number is irrational, and the motion in 
phase space is quasi-periodic. 

When the system Hamiltonian is perturbed such that it becomes non- 
integrable, non-linear resonances appear in the phase space at the location 
of tori with rational winding numbers, altering the topology and destroying 
these tori (see Section H). Trajectories which were confined on the vanished 
tori, now traverse the same general region of phase space. However the 
KAM theorem states that tori with irrational winding numbers (KAM tori) 
are not immediately destroyed by small amounts of non-integrability and 
continue to act as tori in the new phase space. For increasing perturbation, 
the non-linear resonances grow and the KAM tori are eventually destroyed 
by nearby resonances. Within a given region of phase space, those with the 
most irrational winding numbers are the last to break up. Irrationality is 
measured by the rate of convergence of the continued fraction representation, 

w = [ao,ai,a2,...]=ao-\ (2) 
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When a KAM torus is destroyed, a cantorus is left in its place. The 
properties of a broken torus or cantorus in a quantum system have been the 
subject of considerable interest and several numerical studies p3|-[T7|. 

3 Our analytical system 

The atom interacts with a standing wave of near-resonant light (frequency 
ujl) which is temporally modulated with period T. When the detuning 
6l = i^o — (^L (where loq is the resonant frequency of the transition) is suf- 
ficiently large compared to the resonant Rabi frequency 0/2 (proportional 
to the square root of the standing wave intensity), the amplitudes of the 
excited states can be adiabatically eliminated. The dynamics are governed 
by stimulated two-photon scattering between ground states, with momen- 
tum changes in units of 2hkL. Classically the atom behaves as a dipole in 



a conservative one dimensional potential. The Hamiltonian in this limit is 
given by 

n = — — : COS zklx > f { — — n \ (3) 

n=—oo ^ ' 

where f{t/T) specifies the temporal shape of the "kicks," (with < f{t/T) < 
1), ki is the laser wave number, and iloff is the effective Rabi frequency. 
For a two-level atom Q^g = Q^/5l but for our system we instead have 
QcS = ^^(545/(^45 + 844/544 + 843/1543), where the terms in brackets take into 
account the different dipole transitions between the relevant hyperfine levels 
in caesium (F = 4 — > F' = 5,4, 3). In our simulations we assumed equal 
populations of the Zeeman sub-levels, yielding numerical values for the s^j 
of S45 = 11/27, S44 = 7/36, and S43 = 7/108; d^j are the corresponding de- 
tunings. Note that the different magnetic sub-levels will experience different 
AC Stark shifts. For the smallest detuning used in this work this results in 
a 5% spread in the coupling strength. 

We can write the Hamiltonian in dimensionless form as 

2 '^ 

H = ^-kcos(l) y f{T-n) (4) 
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where k = h^cgk'j^T'^ /2M is a measure of the perturbation of the system 
called the kick strength, r = t/T, cp = 2kLX, p = {2kLT/M)px, and H = 
{4:k\T'^ /M)7i. If we consider the dynamics on a single half wavelength of 
the laser beam, the Hamiltonian is that of a driven rotor, with dimensionless 
parameters I = \ and uJq = k, where ujq is the small amplitude oscillation 
frequency of the rotor. In the quantized model, and p are conjugate 
variables with a commutation relation [p, cj)] = —ik, where k = AhkfT/M is 
our scaled Planck constant. 

Experimentally, we use a double pulse kick (see Figure |l| and Section Q). 
For this /(r) the Hamiltonian can alternatively be written as 

H = — — k y^ am cos((/) — 27rmr) (5) 

where am = j^sinc ^ cos ^. (With the sine function defined as sinc(x) = 
sin(x)/x). The rotor is driven by traveling cosine waves, the speed of the 
m^^ cosine wave corresponding to a dimensionless momentum of p = 27rm. 
There will be a primary resonance in the phase space wherever the speed of 



rotation of the rotor matches the speed of a cosine wave |18| . In the reference 
frame in which the cosine wave is at rest, the rotor will be trapped in a 
pendulum potential. This causes a change in the topology of the classical 
phase space at p = pm, for am / 0. For m such that a^ = there is 
a missing primary resonance. The widths of the primary resonance zones 
are given by 6pm = ^V^'mk- The Chirikov condition for overlap is that 
the spacing between the resonances be equal to the sum of the half-widths 
or \pm — Pn\ = 2y/amk + 2yjank. Overlap of adjacent resonances is the 
mechanism for the destruction of KAM tori. We expect tori in the vicinity of 
missing resonances to survive higher kick strengths than all others. A more 
sophisticated analysis must take into account the appearance of secondary 
resonances arising from interactions between the primary resonances. 



4 Experimental Set-up 

Our experimental setup is very similar to that used in the experiments of 
Ammann et al [19,^ and Moore et al [^. Approximately 10^ caesium 



atoms are trapped and laser cooled in a standard MOT powered by two 
diode lasers operating in the infra-red (852nm). The initial trapped cloud 
has a FWHM of ~ 200/im and a temperature of 10— 15pK . The periodically 
modulated potential (described in Section |^ is provided by a third diode 
laser. The beam from this laser passes through an 80MHz acousto-optic 
modulator (AOM) and a single mode optical fibre which spatially filters the 
light and delivers it to the trap. The beam is then collimated and retrore- 
flected from a mirror on the opposite side of the trap to form a standing 
wave potential across the atomic cloud. The calculated beam waist at the 
cloud is 795pm. This potential is temporally modulated via the rf supply 
to the AOM. Our calculated Rabi frequency at the centre of the trap, for 
maximum optical power is Q/2TT = 310 MHz. A reasonably narrow dis- 
tribution in the kicking strength k is produced by the finite widths of the 
cloud and the beam waist of the kicking potential (RMS spread of 6% and 
kmean ~ O.dAkmax)- In the remainder of this paper, k always refers to kmean- 
In the system Hamiltonian, the pulse train is described by /(r). Any 
pulse of finite length will create KAM tori in the phase space of the system. 
These impenetrable barriers to diffusion were observed in recent quasi 5- 



kicked rotor experiments |1£,20,22| and were discussed in a recent paper 
by Klappauf et al |23|. In these studies the effects of the boundaries were 
avoided by tailoring the pulse length to push them into regions of phase 
space beyond the localization length of dynamical localization. However, 



they present a very interesting area of study in themselves and have been 
investigated both numericahy and experimentally |13-17, 2J, 25|. A pulse 



train consisting of single pulses is not the best system for studying KAM 
tori because the only energetic chaotic sea lies within the tori of lowest 
momentum. Energetic chaotic seas aid the diffusion of particles away from 
the partially permeable barrier and are important experimentally in isolating 
the effect of the boundary in phase space. We use a double pulse per kicking 
cycle where the pulse period is T = 25/xs with pulse width a = 1/20 and 
pulse spacing A = 1/10 (see Figure H). This pulse shape gives energetic 
chaotic seas on both sides of the long-lived KAM tori with the smallest 
absolute momentum (see Figure y). To achieve varying levels of spontaneous 
emission, we varied 8^ = ujq —ujl, the detuning of the kicking potential from 
resonance, while simultaneously altering the beam intensity to maintain a 
constant kicking strength. 

5 Transport through cantori 

Classical particles which are no longer confined to their own tori can diffuse 
rapidly even through a recently broken cantorus. The cantor us no longer 
forms an impenetrable barrier in the classical phase space and trajectories 
from both sides of the original boundary can eventually fill the entire region 
between two adjacent, unbroken tori. For quantum particles however, a very 
different behaviour had been predicted. Various numerical studies of particle 
transport through cantori [p^-p!7| suggested that until the phase space area 
escaping through the cantori per kick cycle was ~ h, the diffusion of quantum 
particles would be restricted to a quantum tunnelling like behaviour. Hence 
the prediction was for a distinct difference between the diffusion rates of 
classical and quantum particles. Figure ^a shows a Poincare section of our 
system for a very low kicking strength. The KAM tori at p = ztlOvr and 
ibSOvr are clearly visible as regions of stability between large chaotic seas. 
Figure ^ shows the Poincare section for k ~ 300, which is comparable to 
the kicking strengths used in our experimental work. The tori at p = ztlOvr 
are no longer visible in the phase space at this resolution; however they 
have been shown to still have a significant effect on the dynamics of both 
the quantum |£J,|2^ and classical systems. 

We prepare our atoms so that they initially lie within the p = ztlOvr 
cantori and monitor their subsequent evolution through the boundary. The 
final momentum distributions of the atomic cloud display distinctive 'shoul- 
ders' at p = zblOvr (as shown in Figure ^) due to the inhibition of diffusion 



introduced by the cantori. This KAM locahzation is distinct from the more 
widely studied 'dynamical localization'. The signature of dynamical localiza- 
tion is an exponential lineshape in momentum space which is markedly dif- 
ferent from the box-like distributions and characteristic shoulders observed 
with KAM boundaries. Also, for this experiment, the localization length of 
the system, Ip ~ 170 is considerably longer than the momentum width of 
the 1"** and 2'^"' KAM boundaries at p = iblOvr and iSOvr and hence KAM 
localization occurs before dynamical localization can have an effect. 

6 Decoherence 

The theory of decoherence provides the most recent and so far perhaps the 
most satisfying of a series of attempts to explain the disparities between 
the predictions of quantum mechanics and the everyday experiences of the 
world we inhabit (see Section ||). The field of quantum chaos provides an 
ideal backdrop for a study of the predicted effects of decoherence. It is now 
widely accepted that sensitive dependence on initial conditions - the hall- 
mark of classical chaos - does not occur in closed quantum systems. This 
raises problems for the quantum-classical correspondence (QCC) principle. 
Quantum mechanics must be able to describe the classical limit of chaotic 
behaviour. In this case, employing the limit as ^ — > is not entirely satis- 
factory. Chaotic systems can develop highly complex phase space structures 
in logarithmically short times and hence the small but non-zero value of h 
is an important factor. 

According to the work of Zurek et al [|lO| , [ll|| , these difficulties in restoring 
the classical behaviour can be eliminated by realizing that it is not possible to 
isolate macroscopic quantum systems from their environment. The coupling 
of the extraneous degrees of freedom of a system to the environment destroys 
the quantum coherences on a time scale inversely proportional to the degree 
of coupling. In our experiments and simulations, we introduce coupling 
via spontaneous emission induced by the kicking potential. As the level of 
spontaneous emission increases so does the coupling of the system to the vac- 
uum fluctuations which constitute the environment. The predicted effect of 
this increasing coupling is an increase in the transport of quantum particles 
across the cantori as the quantum diffusion rate tends to its classical limit. 
For the purposes of simulation we note that H = p'^/2 — kcoscj) = -ffught 
while the driving potential is 'switched on', and H = p^/2 = //dark other- 
wise. Classical trajectories involve periods of pendulum motion described by 
Jacobi elliptic functions, alternating with periods of free evolution. Usually 



10^ trajectories are followed during each run. For each pulse, the amplitude 
and phase of the elliptic function is matched to the position and motion of 
each trajectory by inverting the elliptic function numerically. To simulate 
our experiments, we choose a thermal random distribution of initial condi- 
tions. To produce Poincare sections, more uniform conditions are chosen, 
and optimized to reveal the phase space structure. 

The quantum system is represented by a basis of A^ = 128 momentum 
eigenstates |n), where p\n) = nk\n) and n = —64,... ,63. The evolution 
operator for a single kick is 

TT t -l^if darks , .-f^lightx / .-f^darkN 

V = exu(—i — -— — exp(— i — ^ expi—i -—] . . . 

^^ 40fe ^ ^^ 20^ ' ^^ 20^ ' 

I .-f^lights , .17-ffdarkN .„n 

where (m| iJjark \'n) = \n^k'^5m,n, {'m\ i^ught 1^) = \n^k'^5m,n - \k{5m,n+i + 
5m,n-i) and we treat n as periodic. The N x N matrices are exponentiated 
numerically. In the case of evolution without spontaneous emission, the 
evolution is entirely coherent, and is solved by finding the eigenvectors of 
the evolution operator (Floquet method). For incoherent evolution, the 
effects of spontaneous emission are simulated by adding the density matrix 
to two shifted versions of itself, once per kick: 

{m\ V \n) = 2?? ( {m + 1| "P |n + 1) + (m — 1| "P |n — 1) I + (1 — r/) {m\ V \n) 

(7) 

where V is the density operator and r/ is the probability for a particular 
atom to spontaneously emit during one kick cycle. 

A convenient way to visualize the information represented by the density 
matrix is in the form of a Wigner function. For a discrete, truncated basis 
we use the toriodal Wigner function as defined in \2t\. 



w{Xk,Pi,t)= Z^expU— I ( 
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(8) 



Where Pi = {k/2)l and X^ = nk/N. This gives a Wigner function defined 
on a 2A^ x 2N grid. Averaging over cells of four adjacent points we reduce 
the grid to N x N. Examples of results obtained are shown in Figure ^ We 
see that decoherence smooths the Wigner function, removing the rapid oscil- 
lations and negative regions which are characteristic of quantum interference 
phenomena. 



7 Experimental results 

In the absence of decoherence the caesium atoms behave as quantum parti- 
cles. Thus even for cantori where significant diffusion of classical particles 
can occur, our quantum particles should still be strongly contained as can be 
seen from our classical and quantum simulations in Figure ^ where we calcu- 
late the percentage of atoms that cross the p = itlOTr cantori as a function 
of the number of kick cycles. Not until the phase space escaping through 
the cantorus per cycle is ~ ?i, or in our scaled units ~ ^, do the quantum 
particles move to any great extent across the boundary. For our experi- 
mental parameters, the phase space flux through the p = ztlOvr cantori per 
kicking cycle is ~ 4.6fc so the quantum diffusion will be strongly inhibited. 
The predicted effect of increasing the coupling to the environment then is to 
increase the transfer of atoms across the cantorus. As the quantum system 
becomes more and more strongly coupled to its environment, the behaviour 
of the atoms is expected to approach the classical limit of rapid diffusion. 

Our experimental results support this prediction (see Figure |^ . For a 
given kick strength, we tuned closer to resonance in steps from 5 = 2.8GHz. 
A lower limit on the detuning is imposed by the approximation made in our 
numerical calculations that the excited state of the atoms can be adiabati- 
cally eliminated (see Section ^). The upper limit {2.8GHz) is mainly due to 
optical power restrictions from our diode laser - as the detuning increases, 
the intensity must also increase to maintain a constant kicking strength. As 
the detuning decreased, the percentage of particles outside the 1** KAM 
boundary increased and the experimental curve rose towards the classical 
prediction. Our results also show reasonable agreement with our quantum 
mechanical simulations which included the effects of spontaneous emission. 
The main sources of error are the measurement of optical beam power and 
the finite resolution of the CCD camera. 

8 Conclusion 

Using laser-cooled caesium atoms, we have observed the controlled decoher- 
ence of a real quantum system via coupling to the environment. This adds 
to the previous work on decoherence through the atom optics realization 



of the (5-kicked rotor and also the experiments of Haroche et al |27| and 
Wineland et al ||2^ . We have demonstrated that the quantum diffusion rate 
tends towards the classical rate with an increasing degree of decoherence. 
The introduction of decoherence via spontaneous emission increases the rate 
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of transport of atoms across the cantori and alters the characteristic shape 
of the KAM locahzed distribution such that it tends towards the classical 
(uniform) distribution. 

The link between the quantum domain and the familiar classical world 
remains a hotly debated topic. The quantum classical correspondence prin- 
ciple requires that quantum mechanics contains the classical macroscopic 
limit. Taking ^ — > is not a realistic strategy in the laboratory. The phys- 
ical effects of a very small, but definitely non-zero, h manifest themselves 
in chaotic systems through highly complex phase space structures of order 
~ h that develop in logarithmically short times ~ ln(l/?i). It appears that 
only through environment induced decoherence can the quantum — classical 
correspondence principle be justified. This interpretation is authenticated 
via the results presented in this present paper, as well as our group's pre- 
vious study of decoherence in the atomic optics manifestation of the delta 
kicked rotor []l9tP0f|. While we can make our ensemble of Cs atoms behave 
like an ensemble of classical particles, we do not infer that there is any de- 
gree of chaos in the presented quantum system. Everything happening to 
our atoms is intrinsically quantum mechanical, but through the slight phase 
shift that an atom's wavefunction acquires when it recoils from the sponta- 
neous emission absorption — emission cycle the momentum localization is 
destroyed. This process can be observed and understood in the laboratory, 



or through Monte Carlo wavefunction |l£l or density of state calculations. 
We can track our theoretical atom's momentum wavefunction as it absorbs 
and re-emits environmental photons, and through this process localization 
is eliminated and diffusion similar to that predicted classically is observed. 
It is still all quantum mechanics; the physical coupling to the environment 
produced behaviour that only "appears" classical. 
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Figure 1: Double pulse showing definitions of a and A 
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Figure 2: (a) is the Poincare section of the system for a very low kick 
strength. The unbroken KAM tori are clearly visible at p = iblOvr and itSOvr, 
(b) shows the same Poincare section for k ~ 300 where the cantori at /? = 
iblOvr have completely vanished from the phase space 
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Figure 3: Momentum distributions for a kick strength k = 280 as a function 
of the number of kicking cycles r. The experimental data (a) with its prob- 
ability of spontaneous emission per cycle r/ = 1.9%, shows its distinctive 
shoulders at p = ztlOvr. (b) is a quantum simulation for the experimental 
parameters with no spontaneous emission, (c) is a quantum simulation in- 
cluding spontaneous emission rj = 1.9% and (d) shows the corresponding 
classical simulation. 
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Figure 4: Wigner functions for the system after 70 kicks, with k = 280. (a) 
corresponds to a probabihty of spontaneous emission per kick cycle f] = 0%, 
giving pure quantum evolution. We see rapid fluctuations and negative 
regions, indicating quantum interference effects, (b) shows results for rj = 
2%, comparable to experimental value. Quantum effects are reduced leading 
to a more classical distribution. 
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Figure 5: Quantum (solid) and Classical (dashed) simulations for a kicking 
strength of k = 270 showing the percentage of atoms to cross the p = 
zblOvr cantori. This graph clearly shows the inhibition to quantum diffusion 
presented by the cantorus. 
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Figure 6: Percentage of particles to cross the p = ztlO-zr cantori. Experimen- 
tal results for rj = 0.0187 (•) and rj = 0.0503 () all at a kicking strength of 
k = 270. Quantum simulations for the experimental parameters are shown 
as dotted lines. The classical simulation (dashed) and quantum simulation 
with zero spontaneous emission (solid) are shown for comparison 
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